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Abstract 
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bremsstrahlung. We consider several magnetic held configurations 
and propagation scenarios for electrons and positrons. We also pro¬ 
vide an improved energy loss function for electrons and positrons in 
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tions, bremsstrahlung losses, ionization losses and Inverse Compton 
losses with an updated InterStellar Radiation Field. 
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1 Introduction 

Dark Matter (DM) constitutes the largest matter component of the Universe, but its nature 
has so far remained elusive. It is searched for in a number of ways, and most notably via its 
possible electromagnetic emission in the Galaxy. 

In particular, an interesting strategy consists in looking for the emissions produced by the 
interactions of relativistic electrons and positrons, injected by DM annihilations (or decays), 
with the galactic environment. These emissions go under the collective name of ‘secondary 
radiation’ and are essentially of three kinds: (i) radio waves due to the synchrotron radiation 
of the e ± on the galactic magnetic field; (ii) gamma rays due to the bremsstrahlung processes 
on the galactic gas density; (iii) gamma rays due to the Inverse Compton scattering (ICS) 
processes on the interstellar radiation field. 

They have received different attention in the literature. Synchrotron emission has been 
considered since a long time in regions close to the Galactic Center (GC), characterized by 
a large intensity of the magnetic field [1]. Its relevance in wider regions of interest in the 
Galaxy has also been highlighted [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], Bremsstrahlung gamma 
rays have mostly been neglected for what concerns DM studies. Recently, however, their 
importance has been recognized, especially in connection with searches for relatively light 
(10-40 GeV) DM signals from the GC [13, 14, 15, 16, 17, 18]. Finally, ICS gamma rays 
have been identified as an important component of the DM gamma ray spectrum mainly in 
conjunction with the models of leptophilic DM featuring a large annihilation cross section, 
proposed in the wake of the lepton excesses measured by Pamela, Fermi, Hess and most 
recently Ams-02 (see e.g. [19, 20, 21, 22, 23, 24] and many subsequent works). 

All these signatures are potentially very relevant and promising. Their practical use, 
however, depends on a number of different choices, e.g. related to the unknown magnetic 
field configuration, to the unknown propagation parameters of electrons and positrons in the 
Galaxy, to the unknown gas distribution, to the unknown profile of Dark Matter etc. Some 
of these uncertainties are also shared with other possible signals from DM in other Indirect 
Detection channels. In pursuing the goal of identifying Dark Matter or better constraining 


2 








it, it is therefore important to be able to perform multi-messenger analyses in a coherent 
framework and therefore to develop a set of coherent, model independent tools. A step in 
this direction is what we would like to make here. 

More precisely, the purpose of this paper is to provide state-of-the-art tools allowing the 
computation of synchrotron and bremsstrahlung radiation for any (weak-scale) DM model, 
for a set of possible astrophysical configurations that bracket the current sensible ranges of 
the uncertainties. This follows the spirit of previous papers ([25, 26]): the concrete goal is 
to enable the ‘DM model builder’ to readily compute the synchrotron and bremsstrahlung 
phenomenology of her model without having to fiddle with the underlying computations or 
even with the intervening astrophysics, but just by choosing which configurations to adopt 
(and being able to adopt the same choices consistently for other indirect detection channels). 
On the way to achieve that, we have to upgrade some ingredients used to accurately compute 
the population of DM-induced electrons and positrons, namely the energy loss function and 
the e ± halo functions (see below). In the spirit of [25, 26] we always employ semi-analytical 
methods rather than fully numerical ones, in order to better control the relevant physics. 

Armed with the tools just described, we will later play the same game ourselves [27] by 
applying them to refine the constraints from synchrotron radiation derived in [3]. 

The rest of this paper is organized as follows. In Sec. 2 we spell out the astrophysical 
ingredients we need. We discuss in some detail the configurations of the magnetic field 
(relevant for synchrotron emission) while we just recall the main points concerning the other 
ingredients (DM distribution, CR propagation parameters, ISRF, galactic gas maps). In 
Sec. 3 we list several new results. First, in Sec. 3.1 and 3.2 we present our improved energy 
loss function and our improved e ± generalized halo functions. Then, in Sec. 3.3 and 3.4 
we present the two main results and numerical outputs of this work: the generalized halo 
functions for synchrotron radiation and for bremsstrahlung emission. In Sec. 4 we conclude. 

2 Astrophysical configurations 

In this Section we spell out the astrophysical ingredients that we use to compute the propa¬ 
gation of electrons and positrons and ultimately the synchrotron and bremsstrahlung signals. 

While there certainly exist some interdependences between the parameters entering in 
these astrophysical ingredients (for instance, the thickness of the CR diffusive halo -called L 
in the notation below- is related to the vertical extent of the galactic magnetic field, since a 
far reaching magnetic field determines a thick CR confinement box) it is beyond our scope to 
impose such correlations by hand. Our aim is to provide the full range of possible choices and 
it is up to the user to choose sensible combinations. In order to consider inter-dependencies 
self-consistently one is better off choosing a fully numerical approach to the propagation of 
cosmic rays, such as by using Galprop [28] or Dragon [29]. 

2.1 Galactic Magnetic Field 

Our Galaxy has a complicated magnetic field structure, and dedicated efforts by several 
groups have been performed in order to map it: for some recent overviews and sets of ref- 

1 We remind that the equivalent tool for ICS radiation is already provided in [25], and will be updated 
soon in order to include the improvements of the present paper. 
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erences see for instance [31, 32, 33, 34]. We recall here the salient features of the inferred 
magnetic field and then define the simplified functional form that we will adopt. 

The total galactic magnetic field B to t is the sum of a large scale regular B reg and a 
turbulent B tuv b component. These, in turn, can be decomposed in different contributions, 
including disk and halo fields. The regular magnetic field is caused by dynamo effects in 
the galaxy and it can be studied with Faraday rotation measurements of nearby pulsars and 
high latitude radio sources, or with measurements of the polarized synchrotron intensity. On 
the other hand, the turbulent magnetic fields are tangled by turbulent gas flows and can 
be traced looking for their unpolarized synchrotron emission. Recent models of the galactic 
magnetic fields have been proposed e.g. in [35, 36, 37, 38, 33]. 

Rather than the detailed magnetic field geography, the overall intensity is more important 
for our purposes. While we keep in mind that the complicated cartography sketched above can 
have an impact on the determination of the energy losses of electrons and on the synchrotron 
emission from DM, we choose to model the disk field strength by a double exponential in z 
and in r, as proposed e.g. by [40] and [41] for the radial part. Namely, we use 

Btot = B 0 exp ^ (1) 

\ JlD zd J 

where r Q = 8.33 kpc is the location of the Sun. We then adopt several configurations for the 
values for the parameters B 0 , ro and Zd, as shown in the table in fig. 1: 

■ Model 1 (MF1 for “Magnetic Field 1” hereafter) is the configuration used in [25] and 
very similar to the one used in the original GALPROP code (it differs by the normaliza¬ 
tion factor Bo, which has changed a few times in the Galprop literature [40, 42, 43]). 

■ Model 2 (“MF2”) is loosely based on the findings of [35] (and previous [41]). Following 
one of the models in [35] we take a value of 2.1 /jlG for the intensity of the disk regular 
field at solar location (we report it to our value for r Q ); we then add an intensity of 3 
fiG to account for the random component. The resulting field is steeper in r and in z 
than MF1 and reaches slightly higher values at the GC. 

■ Model 3 (“MF3”) is modeled following [44], It is substantially higher at the location of 
the Earth and has larger scale heights both in r and in z, i.e. it extends much farther 
out in both directions. 

2.2 Other astrophysical ingredients 

In this section we recall the other astrophysical ingredients involved in the computations. 
We illustrate most of them in Fig. 1. As a general rule, we want to use state-of-the-art but 
standard ingredients, in order to allow easy comparison with other work. 

o The DM density profile in the Galaxy. We adopt the 6 standard profiles as defined 
in [25] (to which we refer for references and some discussion). They always assume 


4 



Angle from the GC [degrees] 
10 " 30 " 1 ' 5 ' 10 ' 30 ' 1 ° 2 ° 5 " 10 ° 20 ° 45 ° 



Dark Matter halo profiles 


Halo 

a 

r s [kpc] 

p s [GeV/cm 3 ] 

NFW 

— 

24.42 

0.184 

Einasto 

0.17 

28.44 

0.033 

EinastoB 

0.11 

35.24 

0.021 

Isothermal 

— 

4.38 

1.387 

Burkert 

— 

12.67 

0.712 

Moore 

- 

30.28 

0.105 


Propagation parameters for e ± 


Model 

6 

1C o [kpc 2 /Myr] 

L [kpc] 

Min 

0.55 

0.00595 

1 

Med 

0.70 

0.0112 

4 

Max 

0.46 

0.0765 

15 


Magnetic field configurations 
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CD 

[kpc] 

ZD 

[kpc] 

MF1 

[40] 

4.78 

10 

2 

MF2 

[35, 41] 

5.1 

8.5 

1 

MF3 

[44] 

9.5 

30 

4 



o 

a. 



z[kpc] 



Figure 1: Collection of the astrophysical ingredients we use. Top row: DM profiles (figure 
taken from [25]) and propagation parameters for electrons and positrons in the Galaxy. Second row: 
magnetic field configurations. Third row: illustration of the ISRF in two sample locations. Bottom 
row: illustration of the galactic gas densities (figure from [13]). 
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spherical symmetry and their functional forms as functions of the galactocentric coor¬ 
dinate 2 r read: 


NFW : 

PNFw( r 

Einasto : 

PEin(r 

Isothermal : 

Piso(r 

Burkert : 

PBur(r 

Moore : 

PMoo( r 



1 + (r /r s f 


Ps 


(1 + r/r s )(l + (r/r s ) 2 ) 



( 2 ) 


The profiles are plotted in the top left panel of Fig. 1 and their parameters reported in 
the corresponding Table. We remind here that they are normalized by requiring that 
the density at the location of the Sun r 0 = 8.33 kpc be 0.3 GeV/cm 3 * * and the total 
mass of the Milky Way be 4.7 xlO 11 M 0 . Satisfying these two criteria allows to fix the 
r s and p s parameters. Other normalizations (e.g. of the density at the Sun) are used 
in the literature and they would modify the profiles. 


o Electron and positron propagation parameters. These have to be plugged in the 
diffusion-loss equation that we will discuss in detail in Sec. 3.2. We adopt the standard 
choices Min, Med, Max as reported in the table in Fig. 1. Here S and /Co are the 
exponent of the energy dependence of the diffusion coefficient and its normalization 
while L is the thickness of the diffusive halo. While an updated assessment of the 
validity of these ranges of parameters would be welcome, especially in the light of the 
wealth of recent data 3 , this would be beyond our current scope. We continue using the 
standard values reported here also for consistency reasons. 

o InterStellar Radiation Field (ISRF). Electrons and positrons propagating in the 
galactic halo lose energy by Inverse Compton scattering on the ambient light. A detailed 
description of this radiation field is therefore important in order to reliably compute 
the energy losses. We adopt the latest radiation maps extracted from Galprop [53]. 
These replace the ones formerly used in the literature, and in particular in [25]. In 
fig. 1 we draw the two maps in two sample locations (at the Earth and near the galactic 
center) and compare them. One clearly sees the three different components (Starlight 
SL, InfraRed IR and the CMB blackbody spectrum). The current map is much more 
detailed and normalization differences of the order of a factor 2 are visible, but the 
overall behavior is confirmed. We will see that these small differences have an (equally 
small) impact on the observables entering e ± propagation. 

o Gas maps. Electrons and positrons also lose energy by processes occurring on the in¬ 
terstellar atomic and molecular gas (Coulomb interactions, ionization, bremsstrahlung). 

2 We distinguish the notation between the galactocentric coordinate r and the cylindrical coordinates (r, z) 
that we will use in most of the following. Obviously r = y/r 2 + z' 2 . 

3 For instance, a string of recent papers, based precisely on synchrotron radio emission [45, 46, 47, 48] but 

also on positrons [49, 50] and somewhat also on gamma rays [51] and antiprotons [52], finds that the thin 

halo predicted by Min is seriously disfavored. 
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We use the gas maps described in [54] and already used in [13]. We refer to the latter 
for some discussion. They are illustrated in fig. 1. The relevant species are atomic (HI) 
and molecular (H 2 ) neutral hydrogen, ionized hydrogen (HII), neutral atomic helium 
(He) and ionized helium (which is however irrelevant for all practical purposes). As 
discussed in particular in [13], these maps represent a reliable description of the coarse 
grained distribution of gas in the Galaxy, but miss important features at small scales. 
In particular, they do not take into account the regions characterized by a much higher 
gas density (up to 2 or 3 orders of magnitude with respect to the coarse grained maps) 
which are known to exist close to the galactic center (typically at r < 200 pc scales). 
For the purpose of the general tools that we are developing in this work, we do not 
correct by hand the coarse-grid maps by adding the high density regions (contrary to 
what was done in [13]) but we will allow the user to change the overall normalization 
of the gas density in the energy loss function that we will describe below. 


3 Results 

3.1 An improved energy loss function for e ± in the Galaxy 

Using the ingredients described above, we compute an improved function describing the 
energy losses of electrons and positrons during their propagation in the Galaxy. It includes 
energy losses by Coulomb interactions with the interstellar gas, by ionization of the same gas, 
by bremsstrahlung on the same gas, by ICS (using the updated ISRF presented in Sec. 2.2) 
and by synchrotron emission, with the choice of the three magnetic held models discussed in 
Sec. 2.1. Schematically: 

dE 

^tot(77, r, z') = ^Coul+ioniz T U>rem T UcS T Ayn (3) 

where E is the energy of the electron or positron and r and z are cylindrical galactic coordi¬ 
nates. Such a function is provided on the website [55] in the format btot [E , r, z , gasnorm, MF], 
where gasnorm allows to change the overall normalization of the gas densities and MF is a hag 
selecting the magnetic held model. We now recall the dif ferent components of this function 
and illustrate its main features in hg. 2 and 3. Details can be found in standard references 
such as [56, 57] as well as in [25, 24], 

• Energy losses by Coulomb interaction and ionization on neutral matter are 
described by 

b neut (E, r,z) = ^ c a T m e Y^ ™iZi ("log E + log (4) 

i \ e 1 / 

where c is the speed of light, ot = 87 rrg/ 3 , with r e = a em /rn ei is the Thompson cross 
section, n* is the number density of gas species i with atomic number Z % and A Ei is its 
average excitation energy (it equals 15 eV for hydrogen and 41.5 eV for helium). 

On ionized matter, one has 

b lon (E, r, z) = ^ c cr T m e n e ("log — + 2 log ( 5 ) 

4 \ m e A pla / 

4 Notice that we will always limit ourselves to the case of relativistic electrons. 
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where n e is the electron density and E p i a = a/47t n e rf m e /a corresponds to the charac¬ 
teristic energy of the plasma. 

The total energy losses for Coulomb interactions and ionization processes, fecoui+ioniz = 
^neut _|_ jjion^ w qj therefore be given by the sum of eq. (4) and eq. (5) with, respectively, the 
densities of ionized and neutral gas species. In both cases, energy losses are essentially 
independent of E , since the constant terms in the brackets are numerically dominant. 


• Energy losses by bremsstrahlung are described by 


^brem(^) f i 


= C 


r 1 

V]ni(r,z) / 

Jo 


d E„E, d(Ji 


7 7 d V 


( 6 ) 


where E 1 corresponds to the energy of the gamma ray emitted in each bremsstrahlung 
process. The differential cross-section reads 


dcri(E,E 7 ) _ 3 a era (JT 
d E-y 8 tt E 1 


1 + 




(7) 


where (j)\ 2 are scattering functions dependent on the properties of the scattering system. 
For a completely ionized gas plasma with charge Z one has 


(£, £ 7 ) = 02 on (£, E y ) = 4(Z 2 + Z) log 


2 E f E — E~ 


m e c* \ E 7 

and thus the energy losses in this regime (‘weak shielding’) read 


Kim = ^flemdTn, Z(Z + 1) flog ^2 ^ ) E. 




( 8 ) 

(9) 


On the other hand, for atomic neutral matter the scattering functions have a more 
complicated dependence, which is usually parameterized in terms of the quantity A = 
kt E g^ e _ E y For the relativistic regime we are interested in, since E > 1 MeV always, 
one basically cares for the limit A —> 0 for which these functions are constant and take 
the following numerical values: 


0?(A = 0) = d” ss = 45.79, 

^2 (A = 0) = 0« ss = 44.46, 

< e (A = 0) = = 134.60, (10) 

02 e (A = 0) = 0g s = 131.40, 

0(i^ 2) (A = 0) ~ 2 0^ 2) ss . 

The subscript ss in this notation refers to the fact that this regime is usually called 
‘strong-shielding’ because the atomic nucleus is screened by the bound electrons and 
the impinging e ± have to force the shield. In this limit the energy losses read 

Cem = ^Ca em a T n; ^0i iSS - ^02 >S sl E ■ 


(11) 













The total energy losses for bremsstrahlung will therefore be given by the sum of eq. (9) 
and eq. (11) with, respectively, the densities of ionized and neutral gas species. In both 
cases, at leading order, energy losses are linearly dependent on E. A further logarithmic 
dependence arises for scattering in ionized medium, while a small additional energy 
dependence is also found in neutral medium if one accounts for the effect of finite A. 

• Energy losses by Inverse Compton Scattering are described, in exact form, by 


fries = 


3c or 




dq n(e) 


(4 7 2 -T £ )g-l 
(1 + T e q) 3 


2q In q + q + 1 — 2 q 2 + 


1 (r.g ) 2 

21 + T e q 


( 1 - 9 ) 




( 12 ) 


where n(e, r, z) is the number density (per unit volume and unit energy) of photons 
of the ISRF, with energy e, 7 = E/m e is the relativistic factor of the electrons and 
positrons and T e = 4e 7 /m e . 

In the Thomson limit (valid for low electron energies), they reduce to the particularly 
compact expression 


fries 


4 C °+ E 2 
3 ml 


de e n(e, r, z ) 


[Thomson limit], 


(13) 


which makes the energy density in the photon bath Uisrf = f deen(e,r, z) apparent. 

The ICS energy losses are proportional to E 2 (as evident in the Thomson expression, 
but also in eq. ( 12 ) noting that 4 7 2 g is the dominant piece at the numerator) for small 
E. For large E, the dependence softens. 


• Energy losses by synchrotron emission are described by 


b 


syn 


4ca T£l2 B 2 
3 ml 8 tt 


(14) 


where B is the strength of the magnetic field. This formula is in close analogy to the 
one for ICS losses: the integral term in (13) and the B 2 term in (14) correspond to the 
energy density in the photon bath and in the magnetic field respectively. In particular, 
synchrotron energy losses are also proportional to E 2 . 


In fig. 3, left panel, we plot the different energy losses discussed above, at the location 
of the Earth. The different dependences on the er energy are clearly shown. Hence, the 
dominant process in the different energy ranges are, in order, ionization (including Coulomb), 
bremsstrahlung, ICS and synchrotron. 

In fig. 2 we plot the total energy loss function in several locations in the galactic plane (left 
panel) and at several galactic altitudes at the location of the Earth (right panel). We compare 
it with the previous version of the same function not including the improvements listed at 
the beginning of this section (dashed colored lines). The main modification is apparent at 
low energies and it is due to the inclusion of bremsstrahlung, ionization and Coulomb losses. 
Being related to the presence of gas, it disappears at the locations outside of the galactic 
disk. 
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e* energy E [GeV] 


Figure 2: Energy loss function for electrons and positrons in the Milky Way. Left panel: in 
the galactic disk (z = 0), at several locations along the radial coordinate r. Right panel: above (or 
below) the location of the Earth along the coordinate z. Here the magnetic field model MF1 has been 
fixed for definiteness. The circled dot identifies the constant value sometimes adopted. The dotted 
colored lines are the same function before the improvements listed in Sec. 3.1. This figure replaces 
the analogous one (fig. 5) of [25]. 


The modifications due to the use of the new ISRF is minimal and mostly concentrated 
at low energies, so it is hidden by the dominant bremsstrahlung, ionization and Coulomb 
losses in most cases except well outside of the plane where the absence of gas makes it indeed 
visible (see the slight difference between the solid and dashed purple lines corresponding to 
z = 15 kpc in the right panel). 

While in fig. 3 left and in fig. 2 we have chosen the MF1 for definiteness, in fig. 3 right we 
explore the impact of changing the magnetic field model. Not surprisingly, in (r, z) = (3, 0) 
kpc the synchrotron energy losses are larger than at (r, z ) = (8.33,1) kpc, and the ordering 
reflects the intensity of the magnetic field in the corresponding model (see fig. 1). 

In the next subsection we employ this improved energy loss function to compute the halo 
functions for electrons and positrons in the Galaxy. 
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Figure 3: Left panel: the different processes contributing to the energy loss function, at the 
location of the Earth. Right panel: the dependence of the energy loss coefficient function on the 

choice of magnetic field model, in two locations. 


3.2 Revised halo functions for e ± in the Galaxy 

We recall that the number density f(E,r,z) of electrons or positrons at the position (r,z) 
per unit energy E is obtained solving the standard diffusion-loss differential equation 

( t —) \ S r\ 

GeVj V ‘f-gs(b(E,r,z)!)=Q(E,r,z), (15) 


where the first term (which accounts for diffusion) is expressed in terms of the propagation 
parameters discussed in sec. 2.2. The source term Q reads 


Q = < 



dE 


(annihilation) 


(decay) 


(16) 


The function dN e ±/dE is the electrons or positron spectrum from DM annihilations in a 
given final state channel /. 

The solution for /, or rather for the energy spectrum of electrons or positrons d& e ±/dE, 
can be cast [25] in terms of a convolution of the injection spectrum dN e ±/dE with the 
generalized halo functions I(E, E s ,r, z), which are essentially the Green’s functions from a 
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Figure 4: Generalized halo functions for electrons or positrons, for several different values 
of the injection energy E s (color coded). This figure replaces the analogous one (fig. 6) of [25]. 
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Figure 5: Fluxes of electrons or positrons, after propagation, for the case of annihilations 
(top row) and decay (bottom row), shown at two different locations. 


source energy E s to the energy E : 

dQ P ± , „ . c ^ 

—(S,r, 2 )^-/(S,r, Z ) = 


P 


4n b(E, r, z) 


2 V Mdm 

P 


Mi 


DM 


Er 


m dm 


M dm 


dE s * (E s ) ■ I(E, E s ,r, z ) (annihilation) 
dE 


f 


dE, 


f 


dE 


-(E s ) ■ I(E,E s ,r,z) (decay) 


(17) 
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The halo functions /, which replace those in [25], are available on the website [55], in 
the format ElectronHaloFunctGalaxyAnn[halo,propag,MF] [logioX,logioE s ,r,z] (and 
analogously ElectronHaloFunctGalaxyDec for decay) where x = E/E s . 

These functions, particularized at the location of the Earth, are plotted in fig. 4 for 
reference. Comparing with the equivalent functions presented in [25], the main difference 
consists in the evident rise towards small values of the electron energy fraction x , 5 which 
is the direct consequence of the additional, low-energy losses. For small injection energies 
(warmer colors), the rise occurs ‘early’ while moving towards small x, consistently with the 
fact that the new losses are already relevant. For large injection energies (cold colors) the 
rise occurs at small x when E ~ 10 GeV (the regime at which the new losses set in). For 
Min the rise does not happen for large injection energy, as e ± are not efficiently confined on 
the characteristic scale of the energy losses. At a location closer to the Galactic Center (not 
plotted), where energy losses are more relevant, the rise is present. 

In fig. 5 we show the electron spectra, for a few cases. These are in direct correspondence 
with the left panels of fig. 13 of [25]. The differences amount to a factor of a few, up to 
almost one order of magnitude, especially at low energies (where indeed the new losses are 
effective). In the right panels of fig. 5 we show the spectra computed in a location closer to the 
Galactic Center. It is curious to note that, in this case, the fluxes do not follow the intuitive 
normalization ordering MlN —» Med —> Max; in fact, Max yields the most suppressed flux. 
This is just a consequence of the relative importance of the various propagation parameters 
which is different in the Earth’s local neighborhood with respect to that location. Indeed, 
Min, Med and Max are determined as the sets that minimize/maximize the fluxes at Earth. 

3.3 Synchrotron halo functions 

In this subsection we want to obtain the generalized halo functions for synchrotron emission 
which constitute one of the two main technical outputs of this paper. We first review the 
basics of synchrotron emission and then come to the definition of the functions we need. 

The synchrotron power (in erg s ^ 1 Hz -1 ) emitted in a certain frequency v by an isotropic 
distribution of relativistic electrons with energy E in a uniform magnetic field is 

iP S yn(G E,a) = V 3 ————y—F(x) (18) 

m e c z 

with 

1 3 e 

x = vf z/ , v' c = -v c sin a, v c = --5y 2 . 

2 2tt m e c 

Here B is the strength of the magnetic field, a the angle between the line of sight and the 
magnetic field direction and 7 = E/m e the Lorentz factor of the electron or positron. The 
synchrotron kernel F(x) is 

poo 

E(x) = x K 5/3 (x')dx' 

J X 

where K n is the modified Bessel function of the second kind of order n. In presence of a 
randomly oriented magnetic field, which is the case of our interest, the synchrotron power 

5 The functions are defined down to the value of x corresponding to E = 1 MeV, to avoid the regime of 
highly non-relativistic electrons. 
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has to be averaged over the pitch angle a\ 

i r 

Vsyju, E) = - I da sin(a) V syn (u, E, a) 
For relativistic electrons (7 > 2) this corresponds to [30]: 


e 3 B 

V syn (iy, E ) = 2y / 3 y 2 


rrifC* 


K t / 3 (y)K I/3 (y) - -yfu/ 3 (yf - K l/3 (y)'- 


(19) 


( 20 ) 


with y = u/v c . Integrating this quantity over v yields the total power emitted by an electron 
of energy E in all frequencies, i.e. eq. (14). 

Next, the synchrotron emissivity has to be computed convolving the synchrotron power 
in eq. (20) with the number density of electrons per unit energy f(E,r,z ) (in cm -3 GeV” 1 ) 
discussed in sec. 3.2 


/•M D m(/ 2 ) 

jsyn(G r, z) = 2 / dE P sy >, E) f(E, r, z ) (21) 

J m e 

where the minimal and maximal energies of the emitting electrons are determined by the 
electron mass and the mass of the DM particle. The ‘/2’ notation applies to the decay case. 
The overall factor 2 takes into account that, besides the electrons, an equal population of 
positrons radiates. 

Finally, the observable in which we are interested is the intensity 3 of the synchrotron 
emission (in erg enr 2 s -1 Hz ” 1 sr” 1 ) from a certain direction of observation. This is obtained 
by integrating the emissivity of eq. ( 21 ) along the line-of-sight. Schematically: 

3{v,b,i)= f ds h ^ Z) ( 22 ) 

J l.o.s. 47T 

where it is intended that a point in (r, z) is identified by the the parameter s along the line 
of observation individuated by the galactic latitude b and longitude l\ r(s,£,b), z(s,£,b) 6 . 

Recollecting eq. (22) and eq. (17), the synchrotron intensity 3 at a given frequency v and 
for given galactic coordinates ( 6 , £) can be cast as: 




1 ( q \ 2 r M dm dN^ 

0 / dE s y^ j (crv)f—-j^(E s ) I syn (E s , u, £, b) (annihilation) 


2 V Mi 


DM 


{ n \ /'■ Mdm / 2 dN 1 

UtH / dE • E V S-3r^ E ‘) v, i, b) (decay) 

\Mdm/ Jm e f 


dNl 


with the generalized synchrotron halo function 7 syn (i/, E s , £, b) defined as 


(23) 


-7syn(T/ s , 17 £, 7) 


ds (p(r, z) 


l.o.s. r o V Pq 


dE W^y mE ‘' r ’ z) ' (24) 


s We remind the explicit relations r(s,£,b) = s sin 6, z(s,£, b) = \JrQ + s 2 — 2 r Q cos b cos L 
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Figure 6: Generalized synchrotron halo functions, for the DM annihilation case. The upper 
9 panels correspond to an NFW profile, the lower 9 to Burkert; the columns correspond to a fixed 
propagation model, the rows to a fixed magnetic field model. 
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Figure 7: Generalized synchrotron halo functions, for the DM decay case. 
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Figure 8: Synchrotron signal (temperature) at 45 MHz, plotted against the galactic latitude, for 
several choices of DM profile and propagation scheme. 


where g = 1,2 for the decay or annihilation cases respectively and again implicitly r(s,£,b), 
z(s,£,b). The units of J syn are erg/Hz. 

The synchrotron halo functions J syn are available on the website [55], in the format 
ISynAnnI [halo,propag,MF] [log 10 E S) log 10 z/,£,6] (and analogously ISynDecI for decay). 
They are also plotted, for reference, in fig. 6 for the annihilation case and in fig. 7 for the 
decay case. 

The last step needed in order to make contact with actual radio surveys consists in 
expressing the synchrotron signal in terms of brightness temperature T(y) (in K) which is 
defined as: 



(25) 


with ks the Boltzmann constant. In fig. 8 we plot such quantity for a few different choices 
of profiles and propagation parameters. Although a comparison with previous results (e.g. 
in [3] and [8]) is not possible in full details, we have checked that, removing our additional 
refinements, we recover those previous results in most cases. ' 

Before moving on, we would like to point out that our tools can be adapted for usage 
in a more general way. Notably, if a user is interested in the synchrotron signatures from 
a custom galactic magnetic held configuration, she can employ our electron and positron 
galactic density and fold it with the desired MF. This is formally not self-consistent (one would 
be computing the synchrotron energy losses with one MF configuration but the synchrotron 
emission with a different one), but it can be acceptable for practical purposes in the conditions 
in which the dominant energy losses are due to other processes like ICS and bremsstrahlung 
(corresponding to high energies or to regions where the magnetic held is not too large, see the 

7 We cannot however fully reproduce the dependence on the choice of profile in [8]: we find that the 
synchrotron signal is independent on the choice of profile at large latitudes (as we expect from the self¬ 
similarity of such profiles at large radii) while their plots show a sizable residual difference. 
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discussion in Sec. 3.1). Technically, one needs to compute the quantity f(E, r, z ) as presented 
in eq. (17) using the electron and positron halo functions provided in Sec. 3.2. The quantity 
/ can then just be plugged in eq. (21) and then eq. (22) to compute the synchrotron emission 
3. The custom configuration of the MF enters in determining the corresponding synchrotron 
power with eq. (20). 


3.4 Bremsstrahlung halo functions 


In this subsection, in turn, we want to obtain the generalized halo functions for bremsstrahlung 
emission, the other main technical output of this paper. The computation follows quite 
closely the one for synchrotron in the previous subsection, using also the formalism for 
bremsstrahlung spelled out in sec. 3.1. We summarize here the main ingredients for com¬ 
pleteness. 

In close analogy with eq. (23), the bremsstrahlung differential flux (in GeV _1 cm -2 s -1 sr -1 ) 
reads: 


r i 


^^brem-y _ 1 ^”© 

dE^dQ ~ 


P& 


2 rMv M 


dN f , 


2 V M dm / 

cMdm/2 


dE s Ibrem(E s , E~ n £, b ) (annihilation) 


P© 


J dE s 


Ml 


DM 


(decay) 


(26) 

where now (in analogy with eq. (24)) the generalized halo function for bremsstrahlung /b re m 
(E s , Ery, £, b ), which has units of GeV, is defined as 


Ibrem(E s , Ery, £,b) — 2 Ery / - 


ds / p[r, z) 


1 . 0 .s. r o V Pe 


v r E <> 


dE 


V hrem (E 1 ,E,r(s,d)) 
b(E,r(s,e)) 


I(E, E s , r,z). (27) 


The bremsstrahlung power consists in 


Pbrem(^ 7 , E,x) = cE^^ n i(?) ^ ^ 

i ^ 


(28) 


where n* are the number densities of the gas species and the bremsstrahlung cross section 
was given in eq. (7). 


The bremsstrahlung halo functions /b re m are again available on the website [55], in the 
format IBremAnnI [halo,propag,MF] [logioE s ,logioE 7 ,GW (and analogously IBremDecI 
for decay). They are plotted in fig. 9 and 10 (annihilation and decay cases), for reference. 


In fig. 11 we plot the resulting bremsstrahlung line-of-sight y-ray fluxes, for a few cases. 
The agreement with previous calculations (notably [13]) has been verified. We also cross 
checked with fully numerical computations done using GaramaSky [45, 58]. While the spectral 
shape is in very good agreement, we find a difference in overall normalization of the fluxes 
along lines of sight passing close to the Galactic Center. This is due to the fact that GammaSky, 
like GalProp, corrects the bremsstrahlung emissivities by adjusting the normalization of the 
gas densities in each galactocentric ring, in particular close to the Galactic Center (see [59], 
and [13] for a short discussion). We decide to instead use consistently the same maps for 
energy losses and bremsstrahlung emission. 
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Figure 9: Generalized bremsstrahlung halo functions, for the DM annihilation case. Anal¬ 
ogously to fig. 6, the upper 9 panels correspond to an NFW profile, the lower 9 to Burkert; the 
columns correspond to a fixed propagation model, the rows to a fixed magnetic field model. 
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Figure 10: 


Generalized bremsstrahlung halo functions, for the DM decay case. 
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Figure 11: Bremsstrahlung 'y-ray fluxes for the case of annihilations (top row) and decay 
(bottom row), shown for two different channels. 


4 Conclusions 

In the quest for the discovery of Dark Matter, it is important to exploit all possible strategies. 
In Indirect searches, in addition, it is important to be able to exploit the multi-messenger 
and multi-wavelength nature of the possible signals. We have here focussed on the secondary 
radiations from electrons and positrons and presented several upgraded and new results. The 
upgradings concern: i) an improved energy loss function (Sec. 3.1) which fully includes low 
energy losses (Coulomb, ionization and bremsstrahlung) and ii) the revised halo functions 
for electrons and positrons (Sec. 3.2). The new results consist in: iii) the synchrotron halo 
functions (Sec. 3.3); iv) the bremsstrahlung halo functions (Sec. 3.4). All the results are 
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provided in numerical form on the PPPC4DMID website [55]. 

These state-of-the-art tools allow to compute the secondary radiation signal (synchrotron, 
bremsstrahlung and Inverse Compton) from any arbitrary DM weak-scale model and will be 
precious and hopefully instrumental in the current era of precision DM indirect searches. 
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